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SUMMARY 

A set of thermovlscoplastlc nonlinear constitutive relationships (TVP-NCR) 
Is presented. This set Is unique and has been developed mainly for application 
to high- temperature metal-matrix composites (HT-MMC) and Is applicable to ther- 
mal and mechanical properties. Formulation of the TVP-NCR Is based at the 
micromechanics level. The TVP-NCR are of simple form and readily Integrated 
Into nonlinear composite structural analysis. Results show that this unique 
set of TVP-NCR Is computationally effective. It provides a direct means for 
predicting complex materials behavior at all levels of the composite simula- 
tion; that Is, from the constituent materials, through the several levels of 
composite mechanics, and up to the global response of complex HT-MMC structural 
components. 


INTRODUCTION 

High temperature metal matrix composites (HT-MMC) are emerging as mate- 
rials with potentially high payoffs In structural applications. Realization 
of these payoffs depends on the parallel and synergistic development of (1) a 
technology base for fabricating HT-MMC structures and components, (2) experi- 
mental techniques for measuring their thermal and mechanical characteristics, 
and (3) computational methodologies for predicting their nonlinear thermovlsco 
plastic (TVP) behavior In complex service environments. In fact, development 
of computational methodologies should precede the other two because the struc- 
tural integrity and durability of HT-MMC can be numerically assessed and the 
potential payoff for the specific application can be closely estimated. In 
this way. It Is possible to minimize the costly and time consuming experimental 
effort that would otherwise be required in the absence of a predictive 
capability. 

Recent research at NASA Lewis Is directed towards the development of a 
computational capability to predict the nonlinear TVP behavior of HT-MMC. This 
capability Is schematically depicted In figure 1. AS can be seen In this 
figure the capability consists of several computational modules encompassing 
the material TVP behavior (bottom), composite mechanics (sides), and the finite 
element analysis of structural components (top). The TVP computational module 
consists of mathematical models which formally and explicitly relate the 
dependence of the constituent material properties oh 
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The notation used In equations (1) and (2) Is as follows: 

Ptf.T denotes the current property of Interest 
P 0 Is the corresponding property at reference conditions 
Is the melting temperature 
T Is the current temperature 

1 0 Is the reference temperature at which P 0 is determined 

Sp Is the fracture stress determined at T 0 conditions 

o 0 the reference stress at which P 0 Is determined 
• 

Sf Is an appropriately selected stress rate, for example, the stress 
rate at which penetration occurs during Impact 

o c Is the stress rate at which P 0 Is determined, and 

o Is the current stress rate 

The exponents n, m, and l are empirical parameters which are determined 

from available experimental data or estimated from the anticipated behavior 
of the particular product term. 

The first term on the right side of equation (1) represents the tempera- 
ture dependence, the second represents the stress dependence, and the third 
represents the rate dependence or the time dependence, In part. The other 
part of the time dependence Is through the direct time Integration as will 
be described later. Equations (1) and (2) describe, then, material behavior 
In the temperature- stress- time space. 

Each term on the right side of equations (1) and (2) describes a mono 
tonic functional dependence of P/P 0 from some Initial property value to 
a terminal or ultimate material state. The specific shape of the function 
depends on the exponent as Is shown in figure 2 for a fixed exponent and In 
figure 3 for a fixed T/Tf ratio. By judicious selection of the exponent 
and the initial and terminal values, a variety of functional dependences can 
be simulated using equations (1) and (2). 


The form of equations (1) and (2) makes It convenient (provides direct 
feedback) to select the various parameters so that the functional dependence 
described Is consistent with the physical considerations. For example. It 
Is well known that the melting temperature (1 H ) Is a fundamental parameter 
In metals and that the mechanical properties are "zero" at Tm. Also, the 
stress at fracture for some reference condition Is readily determined by 
simple experiments. In addition, the ultimate value of the stress rate may 
be determined from high velocity Impact penetration tests. It can be seen 
from figure 2 that the P/Pp Increases/decreases very rapidly as the melt 
Ing temperature (or any terminal value) Is approached. Furthermore, the form 
of the equations make It convenient to evaluate the exponents from available 


data since each term Is "Isolated" from the others, that Is, the other terms 
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can be taken at reference conditions and will be unity. The equations are 
computationally effective since each term requires simple substitution of 
values and one exponentiation. 

The form of the TVP-NCR selected has all the desirable features men- 
tioned In the Introduction. These features can be conveniently sunroarlzed 
Into three groups: (1J physical, (2) fundamental, and (3) computational, as 

follows: 

(1) Physical - The constitutive relationships describe dependence on: tem- 

perature, time, stress, stress rate, and complete property degradation 

as the ultimate value Is approached. 

(2) Fundamental - the constitutive relationships are: 

• generic - they are applicable to all constituent material proper- 
ties (fig. 4) 

• evolutionary - they are easily extended to Include additional 
dependence, for example cyclic (mechanical, thermal) 

• Isomorphic - they have the same form for all the properties 

• unified - they are fully coupled from the Initial to the terminal 
material state 

• universal - they are equally applicable to any three constituents 
(fibers, matrix. Interphase) 

• nondlmenslonal - they are normalizable with respect to reference 
and ultimate values. 

(3) Computational - the constitutive relationships are: 

• computationally effective - they only require simple substitution 
and exponentiation 

• easily integrated Into nonlinear composite mechanics and struc- 
tural analysis codes - they can be fully Integrated using only a 
few programming statements 


APPLICATION 

The 1VP-NCR were Integrated Into a special-purpose computer code for 
structural analysis of turbine blades made from HT-MMC (ref. 1). This code 
Is depicted schematically In figure 1. The TVP-NCR describe the constituent 
material properties In the material space as shown at the bottom of the 
figure. Note that the cumulative time* temperature* and stress at the cur- 
rent state are tracked through the Integrated computational process as shown 
In the figure. Once the current properties for the constituent materials 
have been determined, they are used In the various levels of composite mech- 
anics to generate the quantities required for the global structural analysis 
(ref. 1). 
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The Integrated aralysls generates properties at all levels of the com- 
posite behavior simulation. Behavior of the longitudinal (fiber direction) 
modulus (t-ji) of the constituents and ply are shown graphically In 
figures 7(a) and (b) for cases 1 and 6 throughout the mission duration. The 
modulus decreases Initially during the start-up and climb part of the mls- 
sloh. The modulus levels off during the steady state (cruise) part of the 
mission. Finally, the modulus Increases during the landing and engine cut- 
off part of the mission. The significant point Is that the TVP- NCR appear 
to represent the material behavior as would be Intuitively anticipated for 
this type of flight mission. 

The corresponding behavior for the transverse (perpendicular to fiber) 
modulus (E 22 ) Is shown In figures 8(a) and (b) and for the In-plane shear 
modulus ( G-j 2 ) figures 9(a) and (b). The behavior of these moduli Is 

about the same as that for the longitudinal modulus (En). 

the behavior of the longitudinal thermal expansion coefficient (an) 
throughout the mission Is shown In figures 10(a) and (b) for cases 1 and 6. 
This coefficient Initially Increases rapidly (during climb), levels off dur- 
ing cruise* and slowly decreases to about Its Initial value during landing. 
This type of behavior Is to be expected since the thermal expansion coeffi- 
cients increase with Increasing temperature. Note that even though an 
for the matrix for case 6 is about 30 percent smaller than for case 1, the 
coefficients for the ply and the interphase are about the same. This Illu- 
strates, In part, the restraining Influence of the fibers In the HT-MMC 
behavior and In addition the Importance of having 1VP-NCR defined at the 
micromechanics level. The corresponding behavior for the transverse thermal 
expansion coefficient (a 2 2 ) Is shown In figures 11(a) ind (b). The 
behavior of <*22 is similar to that of an. 

The behavior of the three lowest natural frequencies of the airfoil 
throughout the flight mission Is shown In figures 12(a) and (b) for cases 1 
and 6. Note that each frequency behaves somewhat differently. This Is 
expected since each frequency is Influenced differently by the centrifugal 
force. The second and third frequencies decrease during climb, Increase 
during the early part of cruise, remain constant during the major portion of 
the cruise, and Increase gradually to about their Initial value during land- 
ing and engine cut-off. On the other hand* the first frequency increases 
sharply during climb (due to centrifugal force stiffening), levels off during 
cruise. Increases sMghtly during landing (cooling but speed retained) and 
gradually decreases to approximately its initial value (zero- speed) . The 
coupled behavior of these three frequencies throughout the flight mission 
further demonstrates the computational effectiveness of the 1VP-NCR to 
represent the physics of the HT-MMC from the constituent materials level to 
the component global structural response. 

The longitudinal stress (on) behavior throughout the mission In the 
constituents and in the ply Is shown In figures 13(a) and (b) for cases 1 and 
6. The stress In the fiber Increases very rapidly during climb, decreases 
gradually during cruise, decreases rapidly during landing, and decreases 
gradually to a small residual compressive stress at engine cut off. lhe ply 
stress exhibits the same behavior as the fiber stress but Is much lower In 
magnitude. The stress In the matrix Increases (compressl vely) very rapidly 
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during climb, remains compressive during cruise and decreases gradually to a 
residual tensile value during landing and engine cut-off. 

The corresponding behavior for the transverse stress (o 22 ) Is shown In 
figures 14(a) and (b) and that for the Intralaminar shear stress (o 12 ) In 
figures 15(a) and (b). Note there are three different regions (A, B, and C) 
for the matrix and two different regions (B and C) for the Interphase In 
which o 2 2 and °12 are computed. These regions correspond to the Intra- 
laminar regions shown In figure 4. The Interesting points to note are: (1) 

the matrix in the different regions Is subjected to both tensile and compres- 
sive transverse stresses which can be of substantial magnitude, (2) the 
interphase can be subjected to relatively high transverse tensile stresses 
which may cause Interfaclal damage, (3) the fiber Is subjected to very high 
transverse tensile stresses which could cause fiber splitting, and (4) the 
transverse ply stress Is relatively small compared to the stress distribution 
In the constituents. 

Collectively, the local stress behavior demonstrates the computational 
effectiveness of the 1VP-NCR to predict the Instantaneous behavior of the 
constituents at the micromechanics level as well as at the ply (macromech- 
anics) level. Determination of the stress behavior In the constituents Is 
possible because the 1VP-NCR are referred to the constituent material space 
and the formulation Is based at the composite micromechanics level, lhe 
gradual decrease Indicated for all the stresses during cruise Is caused by 
the material thermoviscoplastic behavior and may be thought of as a form of 
creep. 


POSSIBLE EXIENSlONS/i IM11A1 IONS 

lhe 1VP-NCR can be extended to Include thermal cycle and mechanical 
cycle effects, diffusion, other material degradation effects, as well as 
time directly. Though terms for these factors are easily added since they 
will be of similar form. It Is not necessarily clear which of these will 
contribute to Independent material behavior. 

In the direct time Integration analysis, the effects of temperature are 
directly accounted for. Any ratchettlng, for example, wilt be part of the 
residual displacements. Also, residual thermal stresses and other stresses 
are known and constitute a part of the cumulative stress history. On the 
other hand, vibratory stress effects are not accounted for In the direct 
time Integration of the structural analysis. Though these effects can be 
accounted for through the stress rate, vibratory stress effects may Indeed 
contribute to Independent behavior. Diffusion or any other material degrada 
tlon can be incorporated once the type of degradation has been defined. 

Other extensions will become self evident as HT-MMC start being extensively 
applied In environments where limited o no property data are available. 

Some limitations of the 1VP-NCR described herein are that they: (1) 

must be used at the current Instant of time, (2) must be used with a direct 
time Integration nonlinear composite structural analysis, (3) cannot be 
verified experimentally at all levels of the composite mechanics analysis 
and at the very high temperatures, and (4) do not Incorporate Initial tangent 
unloading or possible shakedown In the classical plasticity sense. Whether 
(3) and (4) are serious limitations Is yet to be determined. At this stage 
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development It Is prudent to say that these TVP-NCR must be used 
judiciously In design studies relying mainly on sensitivity analyses and 
other Judgment factors that are appropriate for the specific case. 

CONCLUSIONS 

/ -ri/p SrU? 1 ? 06 ther ^1scoplast1c nonlinear constitutive relationships 

(TVP-NCR) for high-temperature metal-matrix composites (HT-MMC) has been 

frt d J? d fl S prese ; te J* , 7h1s set of TVP-NCR Is of simple form. Is appll- 
cable to all thermomechanical properties. Is fully coupled. Is readily Inte- 

effertiwp t0 !° n 1 S? mp0S ^ te structura1 analyses, and Is computationally 
effective. Applicability and computational efficiency were demonstrated * 

through an application to a HT-MMC turbine blade structural analysis. Sen- 
sitlvlty analyses Indicated that substantial perturbations In the TVP-NCR 
exponents have rather minimal effect oh the global and local response of the 
structure. These TVP-NCR make It possible to trace the history of Hl-MMC 

mJr^J]n5°!! ,p ; f,en J S K^ 0m fabr ] cat1on throu 9 h service and from the composite 
micromechanics to global composite structural response. The TVP-NCR are 

sultabTe for preliminary designs and parametric studies. They should be 
tally 1 verlfled^as^et 65190 app11cat1ons s1nce the V have not been experlmen- 
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TABLE I. - EXPONENT PERTURBATION EFFECTS ON STRUCTURAL RESPONSE 


Case | Constituent | Property 


Fiber 

Matrix 


Fiber 

Matrix 


Fiber 

Matrix 


P M - Mech. Prop* 

Pj - Thermal Prop. 

R - Reference 
I - Increase 
D - Decrease 
Units: Displ. - In: 


Exponent 


0.3 0.6 

.6 .1 

.8 2.8 
.2 .1 


.3 0.6 


Structural response 


Cruise conditions 


Residual 


0.5 1.0 

1.0 .2 

.8 2.8 
.2 .1 


0.1 0.2 

.2 .05 

1.2 4.0 

.6 .2 


0.5 1.0 

1.0 .2 


0.1 

.05 0.015 


0.1 

.05 0.0J6 


Unt't 

Freq. 

-0.22 : 

3650 

-0.16 

3840 

-0.31 

4230 

-0.20 

3470 

-0.23 

4050 

-0.29 

3660 

-0.17 

4040 


-0.04 4380 


Displ . 


0.00016 


0.00016 


0.00012 


0.00010 -0.04 4370 


0.00020 


0.00018 -0. 


Percent 

freq. 


Cruise/ 

res. 


Unt*t - Deg: Freq.. Hz 
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TABLE II. - EXPONENT PERTURBATION EFFECTS ON PLY STRESSES 


Case 

Constituent 

Property 

Chanqe 

Exponent 

Ply stresses 

Percent 

ao* n 

Cruise conditions 

Residual 

Cruise/ 

res. 

N 

M 

IB 


0*22 

0*12 

o*n 

0*22 

0*12 


Fiber 

■9 

R 

0.3 

0.6 

0.1 








m 


■Jjw 

R 

.6 

a 

.05 

29.2 

8.4 

-13.7 

0.4 



o 


Matrix 


KM 

.8 

2.8 

.1 





■HI 



HI 


MmH 

R 

.2 

a 

.05 





B 

IH 


2 

Fiber 

IttH 

R 

ItFli 

0.6 

0.1 










■cl# 

R 

1 9 

.i 

.05 

27.6 

10.2 

-13.6 

0.5 

-0.4 

B 

-5 


Matrix 

■21 . 



2.0 

.i 










Mm 

HI 

H 

.05 

.05 






H 


3 

Fiber 

Pm 

R 

EH 

0.6 

0.1 


■1 









R 

mm 

a 

.05 

32.1 


-13.7 

0.5 

-0.7 


+10 


Matrix 


D 

1.2 

4.0 

0.1 






1 


■ 


Mm 

I 

0.6 

0.2 

0.05 


H 




H 


Hi 


Pm 

I 

0.1 

0.2 

0.1 











D 

.2 

.05 

.05 

30.8 

7.1 

-13.9 

1.9 

-1.8 


+5 




R 

.8 

2.8 

a 








HBlfi 


Mm 

R 

.2 

a 

.05 






H 


5 

Fiber 

Pm 

D 


1.0 

0.1 











I 

1.0 

.2 

.05 

33.2 

10.4 

-16.1 

5.2 

EEH 

-2.8 

+14 


Matrix 


R 

.8 

2.8 

.1 





■ 





Mm 

R 

.2 

.1 

.05 





H 



6 

Fiber 

■■ 

I 

0.1 

0.2 



m 









D 

.2 

.05 


33.7 

Bfl 

-13.7 

1.4 

-1.4 

-0.2 

+15 


Matrix 


D 

1.2 

4.0 









BB 


Mm 

I 

.6 

.2 









D| 

Fiber 

Hm 

D 

0.5 

1.0 

0.1 










Ktl 

I 

1.0 

.2 

.05 

26.7 

u.i 

-13.7 

0.1 

~ 

-0.2 

-9 


Matrix 


I 

.4 

2.0 

.1 










■£■ 

D 

.i 

.05 

.05 









P M - Mech. Prop. 

Pj - Thermal Prop. 
R - Reference 
I - Increase 
D - Decrease 
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Fig. 1. - Integrated nonlinear composite structural analysis. 
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Fig. 15. - Intralaminar shear stress variation during flight mission, 
Node 6; ply 4. 
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